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Phantom cosmology allows to account for dynamics and matter content of the universe tracing 
back the evolution to the inflationary epoch, considering the transition to the non-phantom stan- 
dard cosmology (radiation/matter dominated eras) and recovering the today observed dark energy 
epoch. We develop the unified phantom cosmology where the same scalar plays the role of early 
time (phantom) inflaton and late-time Dark Energy. The recent transition from decelerating to 
accelerating phase is described too by the same scalar field. The (dark) matter may be embedded 
in this scheme, giving the natural solution of the coincidence problem. It is explained how the 
proposed unified phantom cosmology can be fitted against the observations which opens the way to 
define all the important parameters of the model. 
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1. According to recent astrophysical data the (con- 
stant) effective equation of state (EOS) parameter 
w c ff of dark energy lies in the interval: —1.48 < 
w c ff < —0.72 (see very recent comparison of ob- 
servational data from different sources in Q, and 
see also 0). It is clear that standard A-CDM cos- 
mology is in full agreement with observations. Nev- 
ertheless, it remains the possibility that universe is 
currently in its phantom DE phase (for recent study 
of phantom cosmology, see |j. Lj. |(| and refs therein). 
Despite the fact that it remains unclear how decel- 
erating FRW world transformed to the accelerating 
DE universe, one can try to unify the early time 
(phantom?) inflation with late time acceleration Q . 
In fact, the phantom inflation has been proposed in 
0- The unified inflation/acceleration universe oc- 
curs in some versions of modified gravity Jp as well 
as for complicated EOS of the universe Q (for re- 
cent discussion of similar (phantomic) EOS, see 01 
and time-dependent viscosious EOS[ll|). 

In the present paper we consider unified phantom 
cosmology with the account of dark matter. Due to 
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the presence of scalar dependent function in front 
of kinetic term, the same scalar field may corre- 
spond to the (phantom) inflaton at very early uni- 
verse, quintessence at the intermediate epoch and 
DE phantom at the late universe. The recent tran- 
sition from decelerating phase to the accelerating 
phase is naturally described there too. On the same 
time it is shown that both phantom phases are sta- 
ble against small perturbations, and that coinci- 
dence problem may be naturally solved in our unified 
model. The equivalent description of the same uni- 
fied phenomena via the (multi-valued) EOS is given 
too. In the final section we explain how the proposed 
unified phantom cosmology can be fitted against the 
observations which gives the way to define all the 
important parameters of the model. 
2. Let us start from the following action: 

s = J <rv=£{-Lii- l -u{4>)d^d^ 

-V(<f>)}+S m . (1) 

Here u>(4>) and V(<j)) are functions of the scalar field 
4> and S m is the action for matter field. Without 
matter, such an action has been proposed in [T^ 
for the unification of early-time inflation and late- 
time acceleration in frames of phantom cosmology. 
We now assume the spatially-flat FRW metric ds 2 = 



2 



— dt 2 + a(t) 2 2i=i (dx 1 *) 2 . Let the scalar field <p only 
depends on the time coordinate t. Then the FRW 
equations are given by 



-a =p+p 
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Here p m and p m are the energy density and the pres- 
sure of the matter respectively. The energy density 
p and the pressure p for the scalar field <j> are given 

by 

p = \u(^ 2 + v(4>) , v = - V{4>) . (3) 

Combining Q and J3J, one finds 
2 • 

Uj((j))(j) 2 = -H - (p m +p m ) , 

K z 

V(<P)=^(3H 2 +h)-^^. (4) 

As usually p m and p m satisfy the conservation of the 
energy: 
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As clear from the first equation ©, in case without 
matter (p m — p m — 0), when H is positive, which 
corresponds to the phantom phase, lo should be neg- 
ative, that is, the kinetic term of the scalar field has 
non-canonical sign. On the other hand, when H is 
negative, corresponding to the non-phantom phase, 
lo should be positive and the sign of the kinetic term 
of the scalar field is canonical. If we restrict in 
one of phantom or non-phantom phase, the func- 
tion lo(4>) can be absorbed into the field redefinition 
given by <p = J d4>\J 'lo{4>) in non-phantom phase 

or tp — d<j>yJ—Lo(<j)) in phantom phase. Usually, 
at least locally, one can solve <f> as a function of ip, 
4> = (f)(ip). Then the action Q can be rewritten as 
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V(<fi) \+S, 



(6) 

Here V(ip) = V((f)(ip)). In the sign qp of iJBJ, the 
minus sign corresponds to the non-phantom phase 
and the plus one to the phantom phase. Then both 
of oj(4>) and V{4>) in the action (JTJ do not corre- 
spond to physical degrees of freedom but only one 
combination given by V(tp) has real freedom in each 
of the phantom or non-phantom phase and defines 
the real dynamics of the system. The redefinition, 
however, has a discontinuity between two phases. 
When explicitly keeping u{<p), the two phases are 



smoothly connected with each other (kind of phase 
transitions). Hence, the function lu(4>) gives only re- 
dundant degree of freedom and does not correspond 
to the extra degree of freedom of the system ( in 
the phantom or non-phantom phase). It plays the 
important role just in the point of the transition be- 
tween the phantom phase and non-phantom phase. 
By using the redundancy of u) {(/>), in any physcally 
equivalent model, one may choose, just for example, 
io{4>) as oj{4>) — ujq {<p — <fio) with constants loq and 
(f)Q. If wc further choose ujq to be positive, the region 
given by cj) > <Ao corresponds to the non-phantom 
phase, the region <j> < (f>o to the phantom phase, 
and the point <j> = <f>o to the point of the transition 
between two phases. 

First we consider the case that the parameter w rn 
in the matter EOS is a constant: w m = p m /p m - (I n 
principle, such dark matter may be presented via 
the introduction of one more scalar field). Then by 
using ©, one gets p m = p m0 a~^ 1+Wm \ Here p m o 
is a constant. If oj{4>) and V(cf>) are given by a single 
function g((f>) as 



V{<t>) = -(3 ff '(0) 2 + ff "(<«) 
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with a positive constant go, we find a solution of (0) 
or given by 



t, H = g'(t) , 




Hence, even in the presence of matter, any required 
cosmology defined by H = g'(t) can be realized by 

CD- 

More generally, one may consider the generalized 
EOS like pjj: p m = -p m + F(p m ). Here F(p) is a 
proper function of p m . Using the conservation of the 

energy (J5J) gives a — aoe 3 1 F <p™> . Let us assume 
the above equation can be solved with respect to p m 
as p m = p m (a). Then if we may choose u)((j>) and 
V(4>) by a single function g(4>) as 



>{<£) = --j/W-FU 



a e : 



V{4>) = i (3 5 '(0) 2 + g"{^)) - Pm (ooe»(«) 

+\f (p m (aoe»<«)) , (9) 
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with a positive constant ao , we find a solution of J2J 
or again: 



= i , # = s'(t) (a = a e 9(t) ) 



(10) 



Hence, any cosmology denned by H = g'(t) can be 
realized by ©. 

Since the second FRW equation is given by 



(11) 



by combining the first FRW equation, the EOS pa- 
rameter w e s looks as 



Wctt = - = -1 
P 
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(12) 



Now it is common to beleive that about 5 bil- 
lion years ago, the deceleration of the universe has 
turned to the acceleration. We now show that the 
model describing such a transition could be easily 
constructed in the present formulation. As an ex- 
ample, we consider the model with constant w m (|7|). 
It is also assumed w m > — 1- Choosing g(4>) as 
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(13) 



we obtain 

W(0) = 



4 1 1 

W m + 1 (t S ~ <t>) 2 

o 9o -o , 



v{4>) 
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& [3(w m + If \<j> t s - 
2 ( I 1 
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3 (w m + 1) \ <t> 2 (t s - 
2 




H o — 5o — -. 



Then the Hubble rate H is given by 
2 /l 1 



H 



3 (w m 4- 1) \t t s — t 



Since 



H 



3 (w rn + 1) I i 2 (t s - f) 2 



(14) 



(15) 



(16) 



the EOS parameter w e ff ljT2*|l goes to w m > — 1 when 
t — > and goes to —2 — «; m < — 1 at large times. 
The crossing w c fr = — 1 occurs when H = 0, that is, 
* = i s /2. Note that 



16t, 



27 (w m + l) 3 (t, - t) 2 i 2 
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(3w m + 1) t s 



(17) 



Hence, if w m > —1/3, the deceleration of the 
universe turns to the acceleration at t — t a = 
(3w m + 1) t s /A. The energy density of the scalar 
field 4> and that of the matter are given by 

\ 2 



■it: 



3k 2 (w m + i f (t s - tf t 2 
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(18) 



Then if the coincidence time t c is defined by p\t=t c = 

r 21 ^ 4 

Pm|t=t c , we find i c = t s - V2t 2 /3g K 2 (w m + 1) j 
We may assume i s could be of the order of the age 
of the universe, t s ~ 10 10 yr^ (lO _33 eV) . On the 

other hand k ~ (l0 19 GeV) _1 - (l0 28 eV)~\ Then 
there is a mixing of very large parameter k and small 
one 1/tg in (|14|l , which can be unnatural. 

We now show that the above problem of the un- 
naturalness could be also avoided in the present for- 
mulation. We now propose the second example in 
J7J) given by 



g U) = -alnl-/31n- 

K 



(19) 



Here a and are dimensionless positive constants. 
As we explain soon, we choose ~ O (l0~ 2 ). Note 
that the parameter order of 10 -2 is not unnatural 
since, say 7r 4 ~ O (10 2 ). Ea. (|19|l gives the following 
expressions for u>(<f>) and V(4>): 



2a(3 (p-l + 01n£) 
k 2 (l -/31n^) V 



v(4>) = 



a(3 (3a0 + 0- 1 + /51n^ 



3(lD m + l)Q 



k 2 (1-0 In 2 



, Wm ~ 1 11 n 



3(«j m +l)a 



(20) 
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Supposing go ~ 0(k~ 2 ), there does not appear 
small parameter like l/t s in J3UJ. Now the Hubble 
rate is given by H — a(3/ (1 — /31n (t/n)) t, which is 
positive if 



< t < t. 



Kef 



(21) 



and has a Big Rip type singularity at t — t s (which 
apparently may not occur due to account of quan- 
tum effects 6J). Since 10 61 ~ e 140 , with the choice 



(3 ~ 1/140, we obtain t s ~ KX 10 61 
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-33 , 



■v) , 



whose order is that of the age of the present uni- 
verse. Then now, due to the property of the expo- 
nential function (or logarithmic function), the small 
scale like t s appears rather naturally. We should 
also note that if a/3 ~ O(10 0-2 ) and t is a present 
age of the universe t ~ (l0 _33 eV) , the observed 
value of the Hubble rate H ~ 10~ 33 eV could be also 
reproduced. Since 



a _q/3 2 (ml + a +l-l) 
a (l-/?ln±)V 



(22) 



the universe turns to the acceleration from the de- 
celeration when 



t = t a = Kef a 1 <t s 



(23) 



Since the energy density of the scalar field 4> and 
that of the matter are given by 

3a 2 /3 2 



k 2 (l-/31n£) t 2 



go l-/31n- 



t 



. \ 3(tO m + l)a 



Pm = go 1 - P In 



3(«) m +l)c 



(24) 



the coincidence time t c could be given by solving the 
following equation: 



1-/3 ln^ 

K 



3(tu m +l)a+2 



3a 2 /? 2 
K 2 go 



(25) 



One may regard p m as the sum of the energy density 
of usual matter, like baryons, and that of (cold) dark 
matter. If p corresponds to the energy density of the 
dark energy, the current data indicate that p : p m ~ 
7:3. Then in the present universe, it follows 



1-/3 In- 



t 



3(tu m + l)a+2 



9a 2 /? 2 
10k 2 9o 



(26) 



Thus, in the model (|19|l . the acceleration of the 
present universe and the coincidence problem seem 
to be explained rather naturally. 

Besides the present acceleration of the universe, 
there was a period of the accelerated expansion of 
universe, which is the inflation of the early universe. 
One can further extend the model i|19fl . as the third 
example, to explain on the same time also the in- 
flation of the early universe. By introducing a new 
dimensionless positive constant 7, the following g(<f) 
can be proposed: 



,(,0 = -nl, 1 ( I -fin ( 7 +^J 



(27) 



As in l|19|) . a and (3 are dimensionless positive con- 
stants and it is assumed /?~0(lO~ 2 ). In the limit 
of large t (t 2 / ' k 2 7) being still less than t s in lt?TJl . 
g(t) coincides with that in (|19fl . Since the scale fac- 
tor a is given by a — a^e 9 ^ as in JSJ), the universe 
is invariant under the time reversal t — > — t in l|27|l . 
Then the universe is shrinking when t < and ex- 
panding when t > 0. The scale factor has minimum 
at t — as a = (1 — (/3/2)ln 7 ) a . In the model 
l|27(l , the Hubble rate is given by 

m = , „ , a3 ,.„ N , ,.„ . (28) 



(l-fln(7+^))(7+^) 



Since 



ad 



1-^ln 
2 



0(1 + a)t 2 



7 + 



(29) 



if t > 0, there are two solutions of a = , one 
corresponds to late time and another corresponds 
to early time. The late time solution of a = is 
obtained by neglecting 7 and coincides with 12311 : 
t = ti ~ Ke 1//3 " Q_1 . On the other hand, the 
early time solution could be found by neglecting 
(3, which is O (l0~ 2 ), to be t = t e ~ K-y/j. Then 
the universe undergoes accelerated expansion when 
< t < t e and ti < t < t s . Here t s is Rip time: 
t s = fty— 7 + c 2 /' 3 ~ Ke 1 /' 3 . f e may be identified 
with the time when the inflation ended. One is 
able to define the number of the e-foldings N e as 
N e = In (a (t e ) /a (0)). Then we obtain 



-a In 



l-fln(2 7 A 
.1-|M7) ) 



(30) 
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It is known that N e should be equal or larger 
than 60. We should note that In- function in (|30fl 
cannot be so large naturally since this requires 
(l-(/3/2)ln(2 7 ))/(l-(/3/2)ln( 7 )) ~ e 60 ~ 10 25 . 
Then the ln-function in (|30[1 should be of order of 
unity, which requires that the parameter a should 
be a ~ 10 2 . For example, since we can rewrite (|3U|1 
as 7 = exp(2//3- In 2/(1 - e"^/")) , we find that, 
when a = 1/ (3 = 240, we have N e = 60 if we choose 
7 as 7 = 0.043925 • • • . In the same way one can 
define number of e- foldings in other similar models. 
For instance, in |12|. without matter, the following 
unified model has been considered: 



/(</>) ee g'(ct>) = h 2 



1 



1 



+ ti 



(31) 



Here ho, to, and t\ are positive constants. It is as- 
sumed to > ti . It has been found 0] that H has two 
minima at t = t± = ±\/ (Vq — t\) /2 and at t = 0, H 
has a local maximum. Hence, (early and late) accel- 
erating phantom phase occurs when t_ < t < and 
t > t_|_. The number of e- foldings may be defined as 



N,. : 



In 



o(0) 
a(t_) 




~ I Arctan I t x 



Since In and Arctan-functions should be the order 
of unity, we find ho/t\ (Note that we have assumed 
to > h > 0.) should be O (lO 2 ) so that N e could be 
equal or larger than 60. For example, when to ^S> t\, 



we find N e 



60, we 



Then if we choose 
find N e ~ 60. Hence, it is demonstrated that scalar 
field may play the role of phantom inflaton in the 
early universe and phantom DE in the late universe 
even in the presence of matter. In the intermediate 
phase of the universe evolution the scalar has the 
standard canonical sign for kinetic energy. 
4. It is interesting to investigate the stability of the 
solution (jHJ. It is easier to work without matter, that 
is, we put go in Q or F in © to be equal zero. One 
defines d/dN = (H- 1 ) d/dt, X = fa Y = f(fa)/H. 
Note that A = Y = 1 in the solution (JSJ. By using 
the first FRW equation ©, we find, for the solution 



i-y 2 



H 



(33) 



1 — A 2 

Then by using the FRW equations , © with 



or (JSJ, and the scalar field equation 

= Lu{cj>)4> + \^\fa)^ + mu{fa)4> + V'{fa) , (34) 
one finds 

Consider the perturbations from the solution A = 
Y = 1: X = 1 + 5X, Y = 1 + SY. From §5$, it 
follows 



SX 

dN \ SY 



= M 



SX 
SY 



, M 



-1 1 
-fi -[i 



(36) 

If the real parts of all the eigenvalues of the ma- 
trix M are negative, the solution X = Y = 1 
is stable. The eigenvalues X± are given by A± = 

+ /i) ± + /i) 2 - 8/i } /2. Then the solu- 
tion (JSJ) is stable if and only if [i > 0. From l|33l) . it 
is seen the positive fj, means positive H . Hence, in 
the phantom phase (H > 0), the solution is stable 1 . 

One can check what kind of the EOS for the scalar 
field appears. For simplicity, the universe without 
matter is considered. For the solution JSJ, by using 
the first FRW equation J2J, we find 

±H 2 = ^f(fa) 2 = p , f{fa)=g'{<t>) , (37) 

which may be solved as 4> = f 1 Here 
/ _1 is the inverse of /, that is, if y = f(x), then 
x = f^ 1 (y). By using (JIJ and Q with p rn = p m = 
go = 0, it follows ui(fa) = —rf'ifa) — P + P- Com- 
bining the above equations, the following EOS may 
be obtained: 



(38) 




1 The scalar model of deceleration/acceleration transition 
has been considered in \\A . As also discussed here, the 
solution JHJ with go = is stable in the phantom phase but 
unstable in the non-phantom phase. As one of the eigen- 
values A-|- becomes very large when crossing w = —1, the 
instablity is high there. In order to avoid this problem for 
above toy model, we may consider two scalar model as in 
Appendix^ In case of one scalar model, the instability be- 
comes infinite at the crossing w = —I point, which occurs 
since the coefficient of the kinetic term u(<j>) in vanishes 
at the point. In the two scalar model, one can choose the 
corresponding coefficients do not vanish anywhere. Then 
we may expect that such a divergence of the instability 
would not occur, which can be checked in Appendix 1X1 



G 



Note that / 1 could be multi- valued function in gen- 
eral. For example, in case of l|13(l . EOS is 



P 



-PT 



(w m + l)p 



ts 
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(39) 



Similarly, EOS for other types of scalar couplings 
may be constructed which shows that scalar field 
dynamics may be always mapped into the (compli- 
cated) EOS. 

5. The description of unified phantom dynamics 
in terms of EOS can be fitted against observations 
if one selects suitable sets of data at low redshift 
(z ~ - 1), medium redshift (1 < z < 100) and 
extremely high redshift (100 < z < 1000). Specifi- 
cally, observational evidences point out that the evo- 
lutionary history of the universe comprises two pe- 
riods of accelerated expansion, namely the inflation- 
ary epoch and the present day dark energy domi- 
nated phase with an intermediate decelerated phase 
where a component of cosmic fluid (dark/baryonic 
matter) has given rise to clustered large scale struc- 
ture. As we have seen, a single (effective) fluid 
may indeed be responsible of both periods of ac- 
celerated expansion. At the same time, this fluid 
should be subdominant during the radiation/matter 
dominated epochs to give rise to baryogenesis and 
structure formation. In any case, whatever the fluid 
is, in order to achieve a unified model which could be 
matched with observations, the cosmological energy 
density has to scale as 



,(«) =a^- 3 (i + (i + 

V a ) \ o-de J 



(40) 



with J\f a normalization constant, (77, \) slope pa- 
rameters and ai <C (lde two scaling values of the 
scale factor which lead to early inflationary epoch 
(1) and late dark energy epoch (DE). It is conve- 
nient to rewrite Eg . I|40fl in terms of the redshift z: 



p{z)=N{l + zf 1 + 



1 



zde V 

1 + zjJ l + z ) 

(41) 

To match this formalism with the above results, we 
have to consider z = a^/a — 1, setting, as standard, 
ao = 1 with the subscript denoting quantities eval- 
uated at the present day, i.e. z = 0. Immediately, 

the Hubble rate is H = f-r. Then, taking into 

2+1 ' & 

account Eq.JHJ), all results in terms of time can be 
translated in terms of redshift. However, the phe- 
nomenological parameters r\ and %, assigning the 
slope of p, are derived approximating the function 



g(t), while the rip time t s and the coincidence time 
t c are, for each given g{t)-vao(le\, related respectively 
to zi and zde- In the same way, using ©, the other 
characteristic times as t a or t e can be translated in 
terms of redshift. Keeping in mind these consid- 
erations, the following discussion holds, in principle, 
for any unified phantom cosmology, comprising dark 
matter, like those which we have discussed. Essen- 
tially, g(t) leads the slope of p and the transition be- 
tween the various epochs. From Ea. l40|) . it is easy 
to see that : 



-("+ 3 ) for 



a<a;< clde 



p^a 3 for a/^a<Ca_D_E • 
, p ~ a x ~ 3 for ai <C olde <C a 



Such an energy density scales as dust matter in the 
range ai <C a <C clde- This means that the fluid fol- 
lows matter along a large part of the universe history, 
while it scales differently only during the very begin- 
ning (a <C ai) and the present period (a 3> clde)- 
Moreover, choosing r\ = —3, the fluid energy den- 
sity remains constant for a <ti ai thus behaving as 
the usual cosmological constant A during the early 
epoch of the universe evolution. Finally, the slope 
parameter x determines how the fluid energy den- 
sity scales with a in the present epoch. Considering 
an effective EOS w c ff = p/p, we get, in terms of 
redshift, 



n 



l + z 



x 



3 12 + z + z// ?>\2 + z + z DE 



1 + ZDE -) . (42) 
'J 



It is worth noting that w e s does not depend nei- 
ther on x nor ™ 2fl£ for high values of z, that is 
for early epochs. On the contrary, these two param- 
eters play a key role in determining the behavior 
of the EOS over the redshift range (0, 100) which 
represents most of the history of the universe (in 
terms of time) and the interesting period for struc- 
ture formation. The role of the different quantities 
(77, x, zi , zde) is better understood considering the 
asymptotic limits of the EOS : lim 2 ^± oc w c ff (z) = |, 
where z — > — 00 refers to the asymptotic future. Set- 
ting 7/ = —3, the fluid equation of state asymptoti- 
cally approaches that of the cosmological constant, 
i.e. w\ = — 1. In general, if we impose rj < —1, we 
get a fluid having a negative pressure in the far past 
so that it is able to drive the accelerated expansion 
occurring during the inflationary epoch. 

It is now clear that zj controls the transition to- 
wards the past asymptotic value in the sense that 
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the larger is z with respect to zj, the smaller is the 
difference between w e g(z) and its asymptotic limit 
77/3. This consideration suggests that Zj has to take 
quite high values (indeed, far greater than 10 3 ) since, 
for z S> zi, the universe is in its inflationary phase. 
The present day value of u> e ff is : 



w = 



V 



3(2 + */) 





n- 


H z DE \ 




n- 


- z DE \ 




V2- 


- z DE ) 




V2- 


- z DE ) 



(43) 

where we have used the fact that zj is very large. 
Being zde > 0, in order to have the present day ac- 
celerated expansion, wq should be negative so that 
we get x > 0. Moreover, depending on the values of 
X and zee, wq could also be smaller than w\ so that 
we may recover phantom- like models. The param- 
eter zjj E then regulates the transition to the dark 
energy -like dominated period. It is worth noting 
that, since p scales with a as dust matter for a long 
period, the coincidence problem is naturally solved 
in any theoretical model endowed with these fea- 
tures, as the unified phantom ones presented above. 
Moreover, although wq could be smaller than -1 as 
for a phantom field, the equation of state asymptot- 
ically tends to a value larger than — 1 (provided that 
— 3<?7<— 1) so that the Big Rip is significantly 
delayed or does not occur at all. In order to take into 
account the dust-matter-dominated era, we have to 
ask for an EOS l|42|) of the form w e s (zm) = 0, which 
gives : 



zm 



VdeX - ( 2 + Vde)v ± y/UDEZ(r),x,ZDE,zi) 



2q 



(44) 



with Dde = 1 + zqe and 

Z{v, X, zde, zi) = yDEV 2 + 2 (2-yDE+2z I )rix+yDEx' 2 

(45) 

It is easy to check that, for reasonable values of 77 
and zj, zm is always a complex number and hence 
the equation of state never vanishes and it is always 
negative. Thus, even if its energy density scales as 
that of dark matter over the most of the universe 
life, dark energy cannot play the same role of matter 
since its equation of state is always significantly dif- 
ferent from null. As a result, we have to include also 
the dark matter in the total energy budget. This in- 
trinsic difference in the behavior of dark energy and 
dark matter equations of state is the reason why the 
former, essentially, contributes to the cosmic accel- 
eration while the latter gives rise to the clustered 
large scale structure. This kind of analysis is partic- 
ularly useful in the case of solution l|39|) where the 
dominance of dark matter or dark energy is ruled by 



the parameters t s and w m . 

Summarizing, we have proposed an unifying ap- 
proach to the problem of inflation, dark matter 
and dark energy in the same theoretical framework. 
As we have seen, it is possible to construct exact 
phantom-like cosmological models where all pecu- 
liar eras of cosmic evolution are achieved. Such an 
evolution can be matched with observations using 
the approach outlined in |j where it is considered 
the dimensionless coordinate distance to Gold SNela 
sa mple Il4l and a dataset comprising 20 radio galax- 
ies [la . ll7| . the shift parameter and the bary- 
onic acoustic peak in the LRG correlation function 
pol ] . By this approach, it is possible not only to show 
the viability of the present unifying phantom model, 
but also to constrain its main parameters which are, 
essentially, the times (i.e. the redshift) of transi- 
tion between the various epochs. As an indepen- 
dent cross check, we can also use a recently proposed 
method to constrain the model parameters with the 
lookback time to galaxy clusters and the age of the 
universe thus obtaining consistent estimates for the 
model parameters 21]. In this case we are not re- 
sorting to distance indicators but to cosmic clocks 
which allow to fix with good accuracy the transi- 
tion scales (at least the one from dark matter to 
dark energy dominated eras). Furthermore, the uni- 
fied phantom model can be compared and contrasted 
with the CMBR anisotropy and polarization spec- 
trum, with the data of the matter power spectrum 
and the growth index 0, 0, . These tests make 
it possible to check the model over different redshift 
ranges than the SNela and radio galaxies data, of- 
fering also the possibility to tighten the ranges for 
the different parameters (in particular z\ and zee, 
or ti with i — {e,c, s,a} in the above "time" de- 
scription). As further remark, this unified phantom 
cosmology can be related to fundamental theories 
as discussed in where generalized holographic 
dark energy can be constructed in this scheme. In 
a forthcoming paper, the detailed comparison with 
observations will be developed and discussed. 



APPENDIX A: TWO SCALAR MODEL 



As seen in (|33I36|I . the solution © with go = 
is stable in the phantom phase but unstable in the 
non-phantom phase. As one of the eigenvalues X± 
of the matrix M in l|36|) becomes very large when 
crossing w = — 1, the instablity is very high there. 
In order to avoid this problem, we may consider two 
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scalar model like 

S = J d 4 Xy 
1 



-^v(x)d t ,xd fi X ~ V((f>, x) 



(Al) 



Here r](x) is a function of the scalar field \- The 
FRW equations give 



Then if 



uj((j))^ +r 1 {x)x 2 = *H 



nt,x) = ^(;3H 2 + H 



(A2) 



V(t,t) = -j (3/(t) 2 + /'(*)) 

the explicit solution follows 

<j> = X = t , = /(t) . 



(A3) 



(A4) 

One may choose that u> should be always positive 
and r\ be always negative, for example 

w(0) - -J* - v^+7W} > o , 

V(X) = -4V^+7W<0- (A5) 

Here a is a constant. We now define a new function 

f(cj>, X ) by 

/(</>, x ) = ~( [ d^{4>) + / ^(x)) , (A6) 



which gives f(t,t) = f(t). If V(cf>,x) is given by 
using /O, x) as 



V(<f>,x) 



K 



(90 



9x 



(A7) 



not only the FRW equations but also the scalar field 
equations are also satisfied: 







uj((/))4> + iw'O)0 2 + 3Hui 
, W(0,x) 



??(x)x + 7}V'(x)x 2 + 3Hri(x)x 



dV{ct>,x) 
dx 



(A8) 



In case of one scalar model, the instability be- 
comes infinite at the crossing w = — 1 point, which 
occurs since the coefficient of the kinetic term ui{4>) 
in Q vanishes at the point. In the two scalar model 
in IIAlfl . the coefficients uj(cf>) and r)(<f)) do not van- 
ish anywhere, as in (|A5|) . Then we may expect that 
such a divergence of the instability would not occur. 
We now check this explicitly in the following. 

By introducing the new quantities, Xd,, X x , and 
Y as 



X x ^ x , (A9) 



the FRW and the scalar Eas. (|A8|l are rewritten as 
follows: 



dX 4 
dN 
dX 



v 



dN 
dZ 
dN 



2Hu{4>) 

v'(x) 



(X 2 - 1) - 3(A - Y) , 
(X 2 - 1) - 3(X X - Y) , 



{Xd, (X^Y - 1) 



2H 2 

+X X (X X Y - 1)} 



(A10) 



Here d/dN = H~ 1 d/dt. For the solution Xd, = 
X x = Y = 1. The perturbations are considered as 



Xa, = 1 + SXd 



Then 



X, 



1 + 5X- 



x 



Y = 1 + SY . 
(All) 




3 3 



H' 2 



(A12) 



The eigenvalues of the matrix M are given by solving 
the following eigenvalue equation 



= A + 



x A tt 

H 2 



3k 2 / u'{6) 
A + 



2H 2 
3k 2 
2IP 



A 



Hw{<j)) 

v'ix) 



(A13) 



The eigenvalues are clearly finite. Then even if there 
is an instability, it could be finite. More complicated 
models along this line may be presented as well. 
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